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Techni-dilaton (TD) was proposed long ago in the technicolor (TC) near criticality/conformality. 
To reveal the critical behavior of TD, we explicitly compute the nonperturbative contributions to 
the scale anomaly (9£) and to the techni-gluon condensate (aGp„), which are generated by the 
dynamical mass m of the techni-fermions. Our computation is based on the (improved) ladder 
Schwinger-Dyson equation, with the gauge coupling a replaced by the two-loop running coupling 
a(fi) having the Caswell-Banks-Zaks infrared fixed point a*: a(/i) ~ a — a* for the infrared region 
m < (J, < Atc, where Atc is the intrinsic scale (analogue of Aqcd of QCD) relevant to the 
perturbative scale anomaly. We find that — (#^)/m 4 — > const 7^ and (aG^„)/m 4 — > (a/a C r — 
l) -3 / 2 — > 00 in the criticality limit m/A-rc ~ exp(—ir/(a/a CT — 1) 1//2 ) — ► (a = a* \ a CI , 
or Nf /*■ Nj r ) ("conformal edge"). Our result precisely reproduces the formal identity (8£) = 
(/3(a)/4a 2 )(aG 2 „), where /3(a) = Atc qatc = — (2ci!cr/7r) • (a/ct cr — l) 3 ^ 2 is the nonperturbative beta 
function corresponding to the above essential singularity scaling of m/Axc- Accordingly, the PCDC 
(Partially Conserved Dilatation Current) implies (Mtd / 'm) 2 (_Ftd / m) 2 — —^(O^/m 4 — ¥ const ^ 
at criticality limit, where A'Ito is the mass of TD and Ftu the decay constant of TD. We thus 
conclude that at criticality limit the TD could become a "true (massless) Nambu-Goldstone boson" 
AItd/th —¥ 0, only when m/i^TD — > 0, namely getting decoupled, as was the case of "holographic 
techni-dilaton" of Haba-Matsuzaki- Yamawaki. The decoupled TD can be a candidate of dark matter. 

PACS numbers: 11.25.Hf, 12.60.Nz, 12.60.Rc, 14.80.Va 



I. INTRODUCTION 



The conformal/scale- invariant (walking) technicolor 
(TC) characterized by the large anomalous dimension 
7 m = 1 was first proposed 0, as a solution to the prob- 
lem of the Flavor- Changing Neutral Currents (FCNC) 
in TC, based on the pioneering work by Maskawa and 
Nakajima [|[ who discovered nonzero critical coupling, 
a cl .(=/= 0), for the spontaneous chiral symmetry break- 
ing (SxSB) to take place in the ladder Schwinger-Dyson 
(SD) equation with non-running (conformal) gauge cou- 
pling a{fx) = a > a cl -. 1 Subsequently, similar solution to 
FCNC within the same framework of the ladder SD equa- 
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1 The solution of the FCNC problem by the large anomalous di- 
mension was first considered by B. Holdom [J, based on a pure 
assumption of the existence of ultraviolet (UV) fixed point in TC 
without explicit dynamics and hence without definite prediction 
of the value of the anomalous dimension. 



tion was also considered without usage of the concept of 
the anomalous dimension [f|. See for review 0. 

Due to the (approximate) scale invariance, the theory 
also predicted [l|, [2( a techni-dilaton (TD), a compos- 
ite pseudo-Nambu-Goldstone (NG) boson for the spon- 
taneous (and explicit) breaking of the scale symmetry 
of the TC, as a techni- fermion and anti-tcchni- fermion 
bound state. 

Actually, the mass function of the fermion £(Q) of the 
SxSB solution takes the asymptotic form Q 



E(Q) - m 2 /Q (Q > m) , 
which was interpreted as [l[ 
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FIG. 1: Schematic behavior of the nonperturbative /3(a) given 
in Eq. (QJ. 



with A being the cutoff introduced to the SD equation 
and Cf the quadratic Casimir of the fermion of the fun- 
damental representation of the gauge group. 

Eq.©, often called Miransky scaling, implies [§| that 
the dynamical generation of to by the nonperturbative 
dynamics should lead to the nonperturbative running of 
the coupling a = a(A/m) (— > a CT as A/to — > oo) even 
when it is nonrunning (conformal) in the perturbative 
sense: 



da 2a cr ( a 

/3 NP (a)=A— = — 1 

OA 7r \ a cr 



(4) 



with a„ being interpreted as the UV fixed point. Sec 
Fig- EH- Actually, the mass scale of m has never been 
created from nothing but transferred from the "hidden 
scale" A whose effect persists even when it is removed by 
taking the limit A — > oo while tuning a — > cv cr . 

It was then argued (H that the dynamically gener- 
ated mass m is a renormalization-group (RG) - invariant 
quantity, ^ = 0, and is regarded as generated by the 
dimensional transmutation: 



m = (J, ■ exp 



oc(fi) 



da 



(5) 



due 



d<x(u) 



to (nonperturbative) running of a(fj,), such that 
= /3 NP (a(/j)), which reflects the (nonperturba- 
tive) scale anomaly 

An. («) 



4a 2 
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due to the mass generation of to. Note that the non- 
perturbative beta function ((4]) has a multiple zero, i.e., 
(3 ~ — [a/a cr — l) s with 5 > 1, which is crucial for Eq. 
© with Eq. ((4]) to reproduce the essential singularity 
scaling Eq.© for /_t = A. 

Initially it was assumed 0, that 



(aGt 



0(m 4 ) , 



so that 



(O=/3 NP (c0-O(m 4 ) 



(7) 



(8) 



namely (0£}/m 4 ~ /3 NP (a) — > in the criticality limit 
a — > a cr - From the scale anomaly through Partially 
Conserved Dilatation Currents (PCDC) , the TD mass 
(Mtd) and its decay constant (-Ftd) are related as 



Mf D F| D = -4<0£) = - 



(9) 



which would imply (M td /to) 2 (F td /to) 2 ~ — /3 NP (a) — > 
for the criticality a — > a CT . It was then argued Q 
that Mtd/w could be arbitrarily small by tuning a as 
/3 NP (o;) — >• (a — >• a cr ), namely, TD could become a 
true Nambu-Goldstone (NG) boson, Mtd/to — >• 0, in 
the criticality limit a — > a CT (to/ A — > 0). 

Actually, Eq. ([8]) (and hence Eq.©) turned out to be 
false at least in the ladder approximation: In the critical- 
ity limit a — > a cr (to/ A — > 0) the straightforward ladder 
calculation 0] of (6fi) = 4(6>q), with the vacuum energy 
(8q) evaluated through the Cornwall- Jackiw-Tomboulis 
(CJT) effective potential, yields 



. N f N TC 4 
-4^—, — TO 



(10) 



in obvious contradiction with Eq.©. Accordingly, we 
have 



4a 2 



i 



to 4 /3 NP (a) m 4 ^ N p(a) 
in contrast to Eq. ©. Eq.© now reads 



(11) 



/ Mtd \ 2 / -Ptd n 



\ to / y to 



-> const 7^ , (12) 



which implies that there is no massless TD in the criti- 
cality limit, Mtd/ui — > const ^ 0, as far as to/Ftd — ► 
const 7^ 0. 2 



with /3 NP (a) given in Eq.©, where (aG^) is the non- 
perturbative contribution to the techni-gluon condensate gauge theories in ladder-type approximation [IoHl3l 



2 There were several other arguments against the massless dilaton 
in the criticality limit of the (nearly) conformal/scale-invariant 
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Recently a possibility was suggested [3] that the TD 
is relatively light compared with other techni-hadrons, 
though not extremely light: The TD mass may be eval- 
uated at a = a C r m t ne limit m/A — > as 



(13) 



through the old calculation [15( of a scalar bound state 
in the gauged NJL model which well describes the 
conformal/scale- invariant gauge dynamics at criticality 
a — > cr cr where the anomalous dimension 7 m = 1 makes 
the induced four-fermion operator marginal with physical 
dimension d = 2(3 — -y m ) = 4 [T^i- Eq.([T3l) is consistent 
with the ladder calculation Eq. ([TU)) (and hence Eq . (TT21 ) . 

Furthermore, in the modern version |16l4l8| of 
conformal/scale- invariant TC based on the Caswcll- 
Banks-Zaks infrared fixed point (CBZ-IRFP), a* = 
a*(Nf, Ntc) [HI, of the two-loop beta function, the cou- 
pling is almost nonrunning a(/_i) ~ a = a* over the wide 
infrared (IR) region fi < Atc below the intrinsic scale 
Atc which is an analogue of Aqcd of QCD (see Eq . ([30l) ) . 
Based on the SD equation in the (improved) ladder ap- 
proximation, with the nonrunning coupling a in the lad- 
der expression simply replaced by the two-loop running 
coupling a(/_t), we have approximately the same result as 
Eq.(TT])-Eq.([3]) with the cutoff A replaced by A T c (to be 
typically identified with the Extended TC scale Aetc); 
and the nonperturbative beta function Eq.(|4]) as well as 
7 m ~ 1 near criticality. 

In this case such a relatively light TD was also sug- 
gested [lil. l20ll from the result of the straightforward cal- 



culation 211 of scalar bound state mass, 



M TD ^ 1.5 m(~ 4FV) ~ V2m 
< M p ,M ai —4.2 771 , 



(14) 



Nf), where N c f r 



m 



and M p ,M ai 



near the criticality a» ~ a cr (Nf ~ 
N"(Ntc) is determined by a* = a, 
are mass of techni-p and techni-ai mesons. The calcula 
tions are based on the SD equation and the homogeneous 
Bcthc-Salpctcr (BS) equation in the improved ladder ap- 
proximation. 3 Although the result Eq. (fT4| is evaluated 
not in the critical limit a = a* — > a cr but slightly away 
from it, the result seems to indicate Mtd/'ti — > const ^ 
in the criticality limit a = a* — > a cr , which is consis- 
tent with the ladder calculation Eq. (fT2")) . Numerically, 
Eq.(HU suggests QUO] 



M TD ~ 500 GeV 



(15) 



3 These results are compared with those in QCD: [22 
{M P /F^)/{M P /F^) QCD ~ 1.3, {MaJF^/iMaJF^QCD z 
0.86, while (M T D/F7 r )/(M scalar /F^) Q cD ^ 0.38. 



in the typical one-family TC model near criticality (with 
N f ~ 4N TC = 8-12). 

More recently, TD mass was estimated by Haba, Mat- 
suzaki and Yamawaki (23j ] in the hard- wall type (bottom- 
up) holographic approach including the techni-gluon con- 
densate (aG 2 ^). It was found that for fixed value of S 
and 7 m , 4 the TD mass is a monotonically decreasing 
function 



TD 



(r -> 00) , 



(16) 



of the techni-gluon condensate T (normalized by the cor- 
responding quantity of QCD), 



1/4 



(17) 



QCD, 



where F n is the decay constant of the (tcchni-) pion of 
order 0(m). It was argued that the limit T — > 00 is 
realized at the criticality /3 NP (a) — > (a = a* — > a cr ), 
as is seen from Eq. (Ti~Tj) . when the value of (aG^) is 
evaluated through Eq.((6]) by assuming the ladder result 
Eq. ([TTJ)) and the nonperturbative beta function Eq.@. 
This would imply the existence of true (massless) NG 
boson at criticality in contrast to Eq.(|13j) and (possibly) 
Eq. ([T4]) . However, from the ladder result Eq. (TT2"1) . Eq. 
implies that 



m 



TD 



F 







TD 



(18) 



in that limit, namely the holographic TD becomes a de- 
coupled TD whose all couplings are characterized by the 
power of (p/Ftu) with the typical momentum p(^ m). 
This is a new feature of the holographic TD. 

The actual phenomenologically interesting situation of 
TC model building is slightly away from the criticality, 
m/A = m/AsTC — 10~ 3 — 10~ 4 =/= 0, 5 in which case 
we have r ~ 7. This implies mass of holographic TD for 
typical conformal/scale-invariant TC model with Nf ~ 
4N TC as: M 



M TD ^ 600 GeV , 

< M p , M ai ~ 3.8 TeV 



(19) 



Note that in the holography the S parameter 1241 and anoma- 
lous dimension f m are not calculable parameters but arbitrary 
adjustable parameters l2- r | . 

In the actual situation of TC, m is the weak scale m ~ TeV 
and A is identified with the typical scale Aetc of the dynamics 
(like the extended TC (ETC)) transmitting the techni-fermion 
mass m to that of the quark/lepton, i.e., A = Aetc ~ 10 3 TeV. 
Thus m/A ~ 10~ 3 — 10~ 4 which corresponds to /3(a) J> 10 -2 
(for N TC = 2-3) from Eq.© and Eq.Q. 
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for the value of S = 0.1 and j m = 1, in rough agreement 
with Eq. (fj~5j) . 6 

Most recently, on the other hand, Appelquist and 
Bai [13] argued, based on the improved ladder SD equa- 
tion with the two-loop running coupling, that there does 
exist a (non-decoupled) massless TD, Mtd/tu — > 0, in 
the conformal/scale-invariant TC in the criticality limit 
(3(a) -> as a = a, -> a CT (Nf — > N"), based on es- 
sentially the same assumption as in Ref. P, 0] , namely 
Eq.©, which is in disagreement with the ladder calcula- 
tion, Eq.|ini), as noted before. (See also Ref. [H].) Note 
that although the beta function in Ref. [27|is somewhat 
different from that in Eq. ((4]) used in Ref. Q , they both 
vanish at the criticality a* a cr (Nf —> Nf 1 ). 

In view of these subtleties in the literature on the crit- 
ical behavior of the TD near the conformal edge associ- 
ated with the CBZ-IRFP, it is very important to settle 
the critical behavior of (aG 2 ^) and (Ofc) in the calculation 
within the same framework as that relevant to the above 
controversy, namely, literally incorporating the perturba- 
tive two-loop running effects as well the nonperturbative 
effects which produce the dynamical mass m. 

In this paper we shall explicitly calculate the nonper- 
turbative contributions to the techni-gluon condensate 
(aG 2 ^) and to the scale anomaly (Ofi) = 4(6>q) arising 
from the fermion mass generation in the TC near con- 
formality /criticality (conformal edge), a* — > a rr (Nt —> 
N"), based on the "improved ladder SD equation" [291 ] . 
Although the improved ladder approximation with the 
two-loop running coupling as well as the ladder ap- 
proximation with nonrunning coupling is not a system- 
atic approximation and hence not very reliable, all the 
above controversy about the techni-dilaton in the lit- 
erature has been confined to this approximation. So 
our aim of this paper is to resolve the confusion within 
this approximation. We first study analytically the solu- 
tion of the improved ladder SD equation, with the two- 
loop running coupling being approximated by a simpli- 
fied ansatz (solution to the "parabolic" beta function 
/3(a) = — boa(a* — a)), 

<V) = rrv^' (20) 

which agrees with the exact two-loop running coupling 
written in terms of the Lambert's W function in the IR 



The value M p , M ai is essentially determined by the value of S: 
Lower S value corresponds to higher M p ,M ai . The calculated 
S value in Ref. [26ll in the same setting as Ref. |21| is higher 
than S = 0.1, which corresponds to the value of M p , M ai in the 
holography close to that of Eq. (114 tt . 



region [i < Atc responsible for the dynamical mass gen- 
eration (sec text). The result will be checked by the 
numerical solution based on the exact two-loop running 
coupling. 

We then calculate the techni-gluon condensate near 
the conformal edge and show explicitly it behaves as 
(G^J/m 4 - (a/a CI - 1)~ 3 / 2 -> oo (a -> a cr ), which 
is a direct evidence against the assumption of Ref. Q (in 
the ladder SD equation with nonrunning coupling) and 
also that of Ref. |27| (in the improved ladder SD equation 
with the two- loop running coupling). Our result directly 
confirms the estimate of the techni-gluon condensate in 
Ref. [HI which indicates divergence of the techni-gluon 
condensate T — > oo at criticality. 

We also find that the numerical calculation of the vac- 
uum energy with the two-loop running coupling agrees 
with the analytical solution (fT0|) with the fixed coupling, 
(#^) ~ — m 4 , again in contrast to the assumption in 
Ref. and Ref. H3- 

On the other hand, the scale anomaly satisfies the for- 
mal relation, (0£) = (3(a) / (Aa) ■ (G 2 ^) . Hence our results 
imply the beta function near criticality: 

4a(<^) fa \ 5 

The result also confirms the assumption made in Ref. [23| 
where the nonperturbative beta function, Eq. (j4]), as well 
as the ladder result of the vacuum energy was used for 
the nonperturbative conformal anomaly to estimate the 
techni-gluon condensate. 

We thus conclude that the nonperturbative beta func- 
tion arising from the nonperturbative effects of the dy- 
namical mass generation in the IR region (fi < Atc) 
is essentially like Eq. ((4]), FigQ] even in the case of the 
two-loop running gauge coupling set in the SD equation. 
It should be considerably changed from the perturbative 
expression near criticality. In Fig. 1141 we schematically 
depict the conjectured behavior of the beta function in- 
cluding both of the perturbative and nonperturbative re- 
gion. 

Our two-loop results combined with the PCDC re- 
lation, i^| D M| D = — 4(0^) ~ to 4 , suggest F^ u /m 2 ■ 
M^Yy/m 2 — > finite at the critical point, which is the same 
as the nonrunning case (fT2]l ; There is no theoretically 
controllable suppression factor for Mtd/to — > 0, as far 
as Ftd/tti is finite. This contradicts the assumptions in 
Ref. @ and Ref. [27j. However, our results cannot ex- 
clude the possibility (the "decoupled TD", Eqs. ((16|) and 
(HI])) that there might exist a very light TD, Af TD - 0, 
if Ftd/w is quite large, as could be the case in the limit 
of the holographic TD [23| . This decoupled TD may be 
dark matter. 
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This paper is organized as follows: In Sec. II, we de- 
scribe the behavior of the beta function in the two-loop 
approximation. We also introduce the parabolic approxi- 
mation in order to solve analytically the improved ladder 
SD equation. In Sec. Ill, we study the analytical solu- 
tion of the SD equation in the parabolic approximation 
and also analyze the numerical solution with the two- 
loop exact gauge coupling. We show that the approxi- 
mation works well. Then we calculate the techni-gluon 
condensate and the vacuum energy. Sec. IV is devoted 
for summary and discussions. 



II. TWO-LOOP p FUNCTION AND 
PARABOLIC APPROXIMATION 




In this section, we study the running effect of the gauge 
coupling constant in the two-loop approximation. It is 
well-known that there appears the CBZ-IRFP [lH, when 
the number of (techni-)fermions is in a certain range, as 
we will show later. If the value of the CBZ-IRFP a* 
slightly exceeds the critical coupling a CI for the SxSB, 
we can apply such gauge theories to the TC with near 
conformality with anomalous dimension 7 m ~ 1 [l6l - [l8j ]. 
We here employ the approach of the (improved) ladder 
SD equation |29||, with the nonrunning coupling in the 
ladder SD equation simply replaced by the running one, 
this time the two-loop running coupling. Although the 
numerical analysis of the (improved) ladder SD equation 
is rather straightforward, it is not so easy to extract nu- 
merically the critical behavior of the solution. We thus 
approximate the two-loop /3 function into a parabolic one 
and will solve analytically the ladder SD equation. In the 
next section, we will demonstrate that this approxima- 
tion works very well 

Let us study the two-loop renormalization group equa- 
tion (RGE) for the gauge coupling constant a [3fJ: 



jti— a = 13(a) = -b Q a 2 - M 3 , (22) 



with 



and 



bi 



12^ 



b = —(UC A -AN f T R ), (23) 

D7T 



17C 2 A -2N f T R (5C A + 3C F ) , (24) 



where Nf represents the number of flavor and the group 
theoretical factors are 



FIG. 2: Behavior of /3(a) in perturbation. The bold solid and 
dashed curves correspond to the two-loop j3 function (1221) and 
the parabolic one (|36[) . respectively. 



for SU(Ntc) gauge theories. 
When b > and b\ < 0, i.e. 



the CBZ-IRFP a* emerges, 

bo 

a* = — r ■ 
-01 

The analytic form of a(fi 2 ) is also known |3l| : 
2 a* 



l + W(z( M 2 ))' 



(26) 



(27) 



(28) 



where W denotes the Lambert function [32| , which is the 
inverse of xe x , and z is defined by 



1 / (i 2 



2 \ boa, 



(29) 



where the intrinsic scale Ato analogous to Aqcd invari- 
ant under RGE is given by [131 : 



Atc = A* • exp 



H ■ exp 




&oa(/i 2 ) bga 



(30) 



C A = N' 



TC, 



Tr — -, c F 



N 2 C - 1 
2^tc : 



(25) 



with the first term in [ ... ] being the usual one- loop con- 
tribution and the second the two-loop one. We can, of 
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course, rescale Axe freely. Here we chose Axe as 

0.78 a*, (31) 



A T c : ol(h = A TC ) = - 



+ W(e- 1 ) 



which reflects the conformal anomaly associated with the 
perturbative running in the UV region fi > Axe domi- 
nated by the one-loop effects, 



perturbative 



g(oQ 

4a 2 



tG 2 



perturbative 



0(A TC ) , 
(32) 

while keeping the (approximate) conformal symmetry 
(via (almost) nonrunning coupling) in the IR region 
ji < Atc so as to be broken only nonperturbatively by 
the dynamical generation of the techni-fcrmion mass m. 
Actually, the UV and IR behaviors of a(/j, 2 ) in Eq. 
are 



1 



and 



a(/x 2 ) 



bn In ■ 



(m 2 > a tc ), 



(33) 



l + e- 



b a* 



( M 2 « A 2 C ), (34) 



respectively. 

Note that in this paper we are not interested in the 
perturbative part of the conformal anomaly in Eq. ([3"2")l 
and will focus on the nonperturbative contributions to the 
conformal anomaly and the relevant techni-gluon con- 
densate associated with the dynamical generation of the 
mass m in the IR dynamics: [H, [23[ 



(^G^t,) perturbative 5 



— {^Lt/full \®u /perturbative 1 



(35) 



the perturbative conformal anomaly 
fj, 1 perturbative = -0(A^ C ) is associated with the 



where 

perturbative running effects of the coupling in the UV 
region (i > Atc- The quantities defined in Eq. ([33]) are 
similar to those discussed in Ref. [23] ■ 

In order to solve analytically the improved ladder SD 
equation, we would need to simplify the expression of 
a. We thus adopt the approximation (|34j) in all region, 
because it enjoys both of desirable natures, the CBZ- 
IRFP (a a* for /i — > 0) and the asymptotic freedom 
(a — > for n — > oo). This approximation corresponds to 
a parabolic /3 function, 



/3(a) = — 6oa(a* — a), 



(36) 



which can be applied from the IR region to the UV re- 
gion. Although the damping of a in the UV region 
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FIG. 3: Behavior of a. The bold solid and dashed curves 
correspond to the behavior of a for the two-loop /? function 
(122l! and the parabolic approximation (|36[) • We took JVtc = 3 
and N f = 11.85, which yields a* = 0.810. The scale «(^ 2 = 
/Li^r) = a a = 7r /4 is given by /x ci ./Atc = 0.00225. Below 
the scale of the dynamical mass m, numerically m/Axc = 
5.88 x 10 -10 obtained by solving the corresponding ladder SD 
equation with the gauge coupling (|28l) . the techni-fermions 
should be decoupled and thereby the running of a is expected 
to be changed. The dots below /i < m corresponds to this 
expectation. The dash-dotted curve below fi < m is for the 
formal solution of the two-loop /3-function. 



is much faster than the two- loop solution (|28[) . (see 
Fig. |31) it turns out that the critical behavior of the dy- 
namical mass is insensitive to the UV behavior of the 
mass function. On the other hand, the linear approx- 
imation P{a) — — 6oa*(a* — a), which yields a(/i 2 ) 

I — e _1 (^ 2 /A TC ) fc ° Q *), is simpler, but it can be ap- 
plied only in a narrower region. 

Schematic behaviors of the two-loop and parabolic f3 
functions are depicted in Fig. O We also show the 
running effects of the gauge coupling a for both cases 
in Fig. [31 where we took iVxc = 3 and Nf = 11.85 
(a* = 0.810). The parabolic approximation is very suc- 
cessful in the IR region fj, < Axc> while the damping 
of a is quicker than the two-loop one in the UV region 
fj, ^> Axe, a s we mentioned above. (See Fig. [3]) Be- 
low the scale of the dynamical mass m discussed in the 
next section, the techni-fermions should be decoupled 
and thereby the running of a is expected to be changed. 
We depict this expectation by the dots below /i < m in 
Fig. [3] 
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III. ANALYSIS OF THE LADDER SD 
EQUATION WITH RUNNING GAUGE 
COUPLING CONSTANTS 

A. The CJT potential and the improved ladder SD 
equation 

The ladder SD approach is a convenient method to an- 
alyze the dynamical generation of the fermion mass and 



its critical behavior. In order to incorporate the running 
effect of the gauge coupling a, a conventional technique, 
so-called the improved ladder approximation, has been 
widely employed p9l |. We can derive the improved lad- 
der SD equation via the CJT potential Vq jt [H| : 



V C jt{B) = - 



N TC N f 



4tt 2 

•A 3 



dxx < — In I 1 
2 



B 2 (x) 



dxx 



dyy 



B{x)B{y) 



(x + B 2 (x))(y + B 2 {y)) 
I 




9(x - y) + ^% 



(37) 



with the (normalized) gauge coupling X(x), 



X(x 



3CVa(^ 2 = x) 



-in 



(38) 



where x and y denote the Euclidean momenta, the full 
fermion propagator inverse is iSj 1 (p) = A(—p 2 )$ — 

B(—p 2 ), and we took the Landau gauge at which 
the fermion wave function renormalization is unity, 
A(-p 2 ) = 1. Sec Fig.H Although the UV cutoff A is not 
needed for case of the two-loop running coupling which 
is asymptotically free in the UV region in contrast to the 
nonrunning case, we have put an artificial A (— >■ oo) in 
Eq.(f3"T|) only for the numerical calculation, which should 
not be confused with A in the nonrunning case used in 
Eqs.Q and ((4]). The variation of Vcjt with respect to 
the fermion mass function B(x) with x = —p yields the 
improved ladder SD equation 



B{x) 



dy 



yB(y) 
V + B 2 (y) 



X(x) 



( X -y)+Me(y- X ) 



(39) 



B. 



Analytic solution for the improved ladder SD 
equation in the parabolic approximation 



CJT 



-i Tr LnS'j 1 



iTrS f S^ + 



FIG. 4: Effective action for the fermion propagator. The 
CJT potential is defined by Tctt = —Vcjt Jdx 4 . Sf and 
So represent the full fermion propagator and the free one, 
respectively. In the last diagram, the solid line with a shaded 
blob and the wavy line represent the full fermion propagator 
Sf and the gauge boson propagator, respectively. 



/3, because of /3 ~ near x ~ and x ~ oo. Under 
these simplifications, we obtain the following differential 
equation and the two BC's: 



dx 



dx 



-B(x)=0, 



and 



(40) 



The integral equation ([5T)|) is equivalent to a set of a 
nonlinear differential equation and boundary conditions 
(BC's). It is, however, difficult to solve analytically the 
nonlinear differential equation in general. We may adopt 
the bifurcation method [Hj], which yields a more handy 
linearized differential equation. We also ignore xdX/dx oc 



(UV-BC) : x— B(x) 
dx 



(IR-BC) : x 2 — B(x 

dx 



+ £(A 2 ) = 0, (41) 



x=A 2 



0, 



(42) 
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where m is defined by the normalization of the mass func- 
tion, 



B(x = m 2 ) = m, 
and the 1RFP A» and the power factor s are 



A* = 



4tt ' 



(43) 



(44) 



Note that in the limit of s — > oo the gauge coupling 
becomes \(x) — \*9(e~ A 2 -,,-, — x), where the step function 
is defined by 0(f) = 1 for t > 0, 6(t) = 1/2 for t = 0, and 
6(t) = for t < 0. 



and 



s — b a* > 0, 



(45) 



respectively. In the parabolic approximation, the nor- 
malized gauge coupling X(x) is 



X(x) 



A* 



1 + (jZ- 



(46) 



We can analytically solve the differential equation (|40[) 
as follows: 



B(aO 



7/7 



= Pi ^ F 



1 — uj — 1 — w w 

1 ) i + -; -x 

2s 1 2s ' s' 

— 1 + W CJ 

~25~' 2s - ' 7 ; ~ ,T * 



(x > m 2 ), 



(47) 



r 



where F(a, j3, 7; z) represents the Gauss's hypergeomet- 
ric function 7 and we introduced 



(48) 



A 2 ; 



and 



1 



A* 

An 



A C r = 



(49) 



The integration constants C\ and d\ are determined 
through the IR-BC and the normalization of B(x). The 
UV-BC gives the scaling relation. 

The normalization B(x = m 2 ) = m yields 



1 — UJ — 1 — UJ UJ 

ciF[ --5—, 5 ,1 + -;-^ 

Is Is s 



diF{ - 



with 



A 2 



On the other hand, the IR-BC gives 



Cl 



ujF - 



1 -uj 
~ 1 2s~' 



A* 



s + uj 
A* 



s — u> 

1+UJ 



2 s 



-1 — a; 

~ 1+ s 

1 + cj 



(51) 



2s 



2s 

1-UJ 

' 2s 



2 + , x 
s 



,2- 



(52) 



where we used Eq. (150)) . 

In the limit of m <C Atc ; we obtain 



7 If we employ the linear approximation, A(x) = A„(l — 
e~ 1 (x/ A^q) b ), the analytical solution is written by the modi- 
fied Bcsscl functions. 



Cl = 



1 + UJ 
2UJ ' 



di = - 



1 -UJ 
2uj ' 



(53) 



which corresponds to the coefficients of the bifurcation 
solution with the fixed gauge coupling X(x) = A„. On the 
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other hand, the UV-BC in the limit of m 2 < A^, c < A 2 
yields 



TC 



r 1 - a r 2 1 



2s 



(1+C,) 2 



r i 



r 2 ( i + 



(54) 

It is noticeable that the dependence of the UV cutoff A 
disappears. Only when uj is pure imaginary, i.e., 



A* > A c 



(55) 



Eq. ([54]) has a relevant solution, 



In 



;27A' 



TC 



WK 2arctancD 
uj Q 



In 



r l-^ r- i+ 



L r r 2 i+i 



2iCj 



-, (56) 



where 



1. 



(57) 



and n = 1, 2, 3, • • • . It is known that the zero node solu- 
tion n = 1 is the true vacuum [33[. Eq. (|56p yields the 
essential singularity scaling relation, 



Atc = A 



tc • exp 



(58) 



similarly to Eq. Q with replacement of a by a* and A 
by A T c ■ 

The behavior of the mass function B{x) in the super- 
critical region A* > A cr is approximately 



.B(a;<A TC ) Vl + w 2 



B(x > A 2 



TO 



sin 



e 2s 



UJ \ Til 

*ln 
2 vm 



■;?) 



f— V 

V TO 2 / 

( — - ) + arctantD 
Vmv 

mAxc 



A(uj) + A{-w) 



(59) 
(60) 



with 



A(uj) 



r h+f r h-i 



e5JA 



TC 



r 2 i 



(61) 



The behaviors of the mass function in the IR and UV re- 
gions correspond to those with the anomalous dimensions 
7 m = 1 and y m = 0, respectively. In particular, owing 
to the quicker damping of a than the logarithm, there is 
no log correction unlike QCD. On the other hand, the IR 
behavior is the same as that for the fixed coupling. 

In passing, the critical number Nf, which corresponds 
to A rr , is 



Nf = 47V T c 



1 



10 5iV 2 c - 3 



Since the power factor s is 

(lliV T c - 2N/) 2 



s = boa* 



-6[17iV 2 c - N f (5N TC + 3C F )] ' 
at the critical point, it reads 



(62) 



(63) 



b a c 



N TC 



18(iV 2 c - 1) 



3^TC 



12 N TC 
~5 57V 2 C - 3 



For Ntc = 3, they are numerically 



417 
"35" 



11.914, 



107 

560 



(64) 



0.19102. (65) 



We can solve numerically the improved ladder SD 
equation Q39p with the normalized gauge coupling 
The computational technique is described in Ref. [3£ 

We depict the analytical and numerical solutions of 
B{x) in Fig. where we took A^c = 3 and N f = 11.63. 
Although we drastically simplified the integral equation 
([3TJ1) into the linearized differential equation (|40[) with the 
two BC's, we find that the approximation works well. 

The scaling relations in the numerical and analytical 
approaches are shown as the dashed and dotted curves 
in Fig. [51 respectively. We confirmed that the numerical 
solution is unchanged for A/Atc = 10 1,2 ' "• 10 . (It is not 
the case for A = Atc, however.) In the figure, we took 
A/ Atc = 10 5 . The shapes of the scaling relation are 
qualitatively similar. We also find that the analytical 
expression (|56[) is close to the numerical solution for the 
two- loop gauge coupling ([28)) . which will be discussed in 
the next subsection. 



C. Numerical solution for the improved ladder SD 
equation with the two-loop gauge coupling 

We studied the parabolic approximation in the analyt- 
ical and numerical ways, so far. 

Let us now solve numerically the improved ladder SD 
equation with the gauge coupling ([28)) expressed through 




FIG. 5: Behaviors of the mass function in the parabolic ap- 
proximation. The solid and dashed curves correspond to the 
numerical solution of the improved ladder SD equation (|39[) 
with the running gauge coupling (|46p and the bifurcation so- 
lution P?l>. respectively. We took N T c = 3 and N f = 11.63, 
which yields a, = 0.900 and A» = 0.287. 



FIG. 6: Scaling relations in several approaches. The solid, 
dashed and dash-dotted curves correspond to the numerical 
solution of the ladder SD equation ((39} for the two-loop /3 
function, the parabolic approximation, and the fixed coupling, 
respectively. The dotted one is for the analytical expression 
(|56[) in the parabolic approximation. We took JVtc = 3 and 
varied continuously Nf. 



the Lambert function. We calculate the Lambert func- 
tion via the Halley's method }32j . 



w j+1 = 



(Wj + 1) 



(nij +2)(uije™i —z) 
2w~+2 



(66) 



The computational technique for solving the improved 
ladder SD equation is described in Rcf. [36| . The re- 
sults arc depicted in Figs.l5HT51 We also show the results 
for the fixed gauge coupling X(x) — A*. In this case, 
we take the UV cutoff A of the SD equation to Atc- 
We confirmed that the consequences of the fixed gauge 
coupling are consistent with those in Ref. [2l|, where 
X(x) = A*#(A TC — x) was essentially used, instead of 
the two- loop one (|28|) . 

We depict the scaling relation in Fig. [5] We confirmed 
that the numerical solution for the two-loop gauge cou- 
pling fl28j) is unchanged for A/A T c = 10 2, " ' 10 . (It is not 
the case for A/Atc = 10 ' 1 , however.) In the figure, we 
took A/Atc = 10 5 . We find that the numerical values of 
the dynamical mass m/Axc for the two- loop gauge cou- 
pling is smaller than those for the fixed coupling. It is 
amazing that the analytic solution for the parabolic ap- 
proximation is quantitatively close to the numerical one 
for the two-loop gauge coupling. 

In Fig. [7J we show the behaviors of the mass functions 
for the two-loop gauge coupling, the parabolic approx- 
imation, and the fixed gauge coupling. In this resolu- 
tion, we cannot distinguish each other. We did not draw 



here the analytical solution (14"T1) for the parabolic approx- 
imation. Although the behavior is close to the numeri- 
cal one, there is a slight deviation between the analyti- 
cal and three numerical solutions. Compare Fig. [5] with 
Fig. [7J An important point is that we normalized the 
mass function by the dynamical mass B(x = m 2 ) = m, 
not by Atc- Note that the dynamical masses for the 
two-loop gauge coupling, the parabolic approximation, 
and the fixed gauge coupling are numerically obtained 
as to/ Atc = 1-08 x 10~ 4 , 0.845 x 10~ 4 , 13.3 x 10~ 4 for 
Ntc = 3 and Nf = 11.63, respectively. If we had normal- 
ized B{x) by AtC; the three behaviors would thus look 
very different. Owing to this universal nature of the di- 
mensionless mass function normalized by the dynamical 
mass, B(x) /m, the normalized physical quantities such as 
the decay constant F^/m, the vacuum energy (6£)/m 4 , 
and the techni-gluon condensate (G 2 lv )/m A , which are 
determined through the mass function, become insensi- 
tive to the approximations of the running gauge coupling 
near the conformal edge, as we will see later. 

How about the relation between the scale a(/j, 2 = 
/i 2 r ) = a cr and the dynamical mass ml By definition, 
we can obtain the scale /i cr by 



W(z CI ) = AX, - 1, 



with 



Z C r 6 



,2 \ °o 



A 2 



(67) 
(68) 




FIG. 7: Behaviors of the mass function in the numerical ap- 
proaches. The bold solid, dashed and dash-dotted curves cor- 
respond to the numerical solutions of the improved ladder SD 
equation for the two-loop f3 function, the parabolic approxi- 
mation, and the fixed coupling, respectively. We normalized 
each scale by each dynamical mass B(x = m 2 ) — m. We 
took JVtc = 3 and Nf = 11.63, which yields a* = 0.900 and 
A* = 0.287. In this resolution, we cannot distinguish each 
other, however. 



For iVrc = 3 and N f = 11.63, the scale cv(// 2 = (J% T ) = 
a cr is numerically obtained as /Uct/Atc — 0.189. As 
we showed previously, the dynamical mass is to/Atc — 
1.08 x 10~ 4 . When we vary the number of flavor to 
Nf = 11.85, they are much more hierarchical, ^ cr /Axc = 
0.00225 and m/A T c - 5.88 x 10~ 10 . In the parabolic ap- 
proximation, we find a more handy formula, 



Mcr — e 2s Atc 



(69) 



This is apparently much larger than the dynamical mass 
m in Eq.flSU), 



"A 



TC ■ 



(70) 



We show a concrete value for JVtc = 3 and Nf = 11.85 
in Fig. O 

Let us calculate the decay constant F v of the techni- 
pion, which is connected with the weak boson mass. We 
assume that a part of the fcrmion flavor Nf couples to 
the electroweak current. In order to estimate the decay 
constant, we employ the Pagels-Stokar formula [37| . 



2 N D N TC 
4ir 2 



[x + B 2 (x)} 2 



x _d_ 

4 dx 



B 2 {x) 



(71) 



where No denotes the number of fermion doublets which 
couple to the electroweak current. The numerical results 



FIG. 8: Decay constant in the numerical approaches. The 
solid, dashed and dash-dotted curves correspond to the nu- 
merical solution of the ladder SD equation (|39|l for the two- 
loop P function, the parabolic approximation, and the fixed 
coupling, respectively. We took A^tc = 3 and Nd = 1. 



for the two-loop gauge coupling, the parabolic approxi- 
mation, and the fixed gauge coupling are shown in Fig. [5] 
We found that the parabolic approximation works well. 
Note that F„/m ~ 0.41 x y/N D N TC /3 near the critical 
coupling, where we took into account the dependence of 
JVtc an d JVd. Thus, when we fix F^ = 246 GcV, we can 
estimate the dynamical mass as m ~ 1 TcV/\/JVd Atc- 
The non-perturbative techni- gluo n condensate defined 
in Ea. ([35]) can be estimated as [35[ 



(G 



2 ) 



2tt 2 



dxx- 



B 2 (x) 
+ B 2 (x) 



(72) 



See also Fig. [9] Note that after substituting the solution 
of the ladder SD equation, the expression of (G 2 ^) ex- 
plicitly depends only on the mass function B(x) and has 
no ex plic it dependence of the running of the gauge cou- 
pling [35j . Hence we expect that the result is not sensitive 
to the details of the running behaviors of the gauge cou- 
pling: The numerical results for the two-loop gauge cou- 
pling, the parabolic approximation, and the fixed gauge 
coupling are shown in Fig. [TO] Wc found that the be- 
havior of (G^„) is not like {G 2 ^) ~ m 4 as assumed in 
Rcfs. 0,113, but (G^/to 4 - 1/w 3 -> oo in the limit 
of A* — > A cr . Our result directly confirms the estimate 
of the techni-gluon condensate made in Ref. [23[ which 
assumed the ladder result for the vacuum energy Eq.(fTTJ|) 
and nonperturbative beta function Eq.@ in the case of 
the nonrunning coupling. We can show this behavior by 
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FIG. 9: Techni-gluon condensate associated with the genera- 
tion of mass m. 
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FIG. 11: Behaviors of the vacuum energy in the numerical 
approaches. The solid, dashed and dash-dotted curves corre- 
spond to the numerical solutions of the ladder SD equation 
(|39|) for the two-loop /? function, the parabolic approximation, 
and the fixed coupling, respectively. The factor r\ is defined 
by r) = N TC Nf/(2TT 2 ). We took 7V T c = 3. 



the vacuum energy, 



FIG. 10: Behaviors of the gluon condensate in the numerical 
approaches. The solid, dashed and dash-dotted curves cor- 
respond to the numerical solution of the ladder SD equation 
(|39|l for the two-loop /3 function, the parabolic approximation, 
and the fixed coupling, respectively. The factor r\ is defined 
by rj = N T cN f /(2-K 2 ). We took N T c = 3. The cutoff is 
A/Atc = 10 5 for the two-loop /3 function and the parabolic 
approximation. For the fixed gauge coupling, the cutoff is 
taken as A = Atc- 



using the approximation (|59)) . i.e. 



/r 2 v „ NrcNf 



A±, 



dxB 2 (x) 



(73) 



and thus 



(C 2 v „ N TC N f l+tf . fA TC 



N TC N f 
2^ 



(74) 

where we used the scaling relation ([55]) . 

Next to the vacuum energy V = (Off) /4 as defined in 
Eq.(f3"5j). Substituting the solution B so \(x) of the ladder 
SD equation ([39]) for the CJT potential ([37]). we obtain 



N TC N f 
8tt 2 



B ao i), 




f A2 




/ dxx 


ln ^1 


Jo 





BUx) 
x + BUx) 



(75) 



(76) 



where we explicitly wrote the subscript for the mass func- 
tion in order to distinguish the vacuum energy from the 
CJT potential itself. (Of course, also in the expressions 
(HH and ([?2"]) . B(x) represents B{x) = B sol (x).) The 
numerical results for the two-loop gauge coupling, the 
parabolic approximation, and the fixed gauge coupling 
are shown in Fig. 1111 It is clear that the vacuum energy 
normalized by m does not vanish, because the numerical 
calculations shown in Fig. [TT] suggests 



(0£) = AV ~ -0.76 rim 4 , with r\ 



N T cN f 
2ir 2 ' 



(77) 



near the critical coupling. This result disagrees with the 
assumption in Refs. [2I [27\. 

The approximate expression (|59|) suggests that our re- 
sults coincide with the estimate of V for the fixed gauge 
coupling 0: 



AV 



AN TC N 



l m 4 = -0.81 rim 4 . 



(78) 
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FIG. 12: Ratio of the vacuum energy and the techni-gluon 
condensate. The solid, dashed and dash-dotted curves corre- 
spond to the numerical solutions of the ladder SD equation 
(|39|l for the two-loop ft function, the parabolic approximation, 
and the fixed coupling, respectively. We took Ntc = 3. The 
dotted line corresponds to the least-squares fitting by £u) 3 /A*. 
Numerically, we obtain £ ~ —0.026. 



FIG. 13: Instantaneous exponent of the nonperturbative beta 
function with respect to Q. The vertical axis at A» = 1/4 = 
0.25 is the criticality (conformal edge). The solid, dashed and 
dash-dotted curves correspond to the numerical solutions of 
the ladder SD equation (|39[) for the two-loop ft function, the 
parabolic approximation, and the fixed coupling, respectively. 



It is to be noted that this value is also close to the numer- 
ical estimate ([77} in our case. In fact, although the CJT 
potential itself explicitly depends on the running of the 
gauge coupling, the vacuum energy has an explicit depen- 
dence only on the mass function B so \{x) after using the 
solution of the SD equation, and hence only depends im- 
plicitly on the running gauge coupling through B ao \(x). 
(Compare Eq. fl37J with Eq. (JTH).) The IR behaviors 
of B so \(x) for the two- loop running and fixed couplings 
coincide each other, as shown in Fig. [7J (Inside of the 
frame of the figure corresponds to the IR region.) Since 
the UV contribution to the vacuum energy is negligible, 
the vacuum energy (|77|) for the two-loop running gauge 
coupling is almost the same as that for the nonrunning 
one. 

From these analytical and numerical analyses, we con- 
clude that in the vacuum energy there is no divergence 
unlike the techni-gluon condensate, (G 2 „)/m 4 — > oo, and 
also the quantity V/ m 4 does not approach to zero in the 
limit of A* — > A cr , within our approximation. 

Since the formal RG analysis yields [l8[ 



(79) 



the ratio 4V '/ (G^) is closely connected with the f3 func- 
tion. We depict it in Fig. [12] By using the least-squares 



method, we numerically obtain 



4a 



W 

(Gm„) 



£— , with £ = -0.026. 
A* 



(80) 



in the form similar to the nonrunning case. In the case 
of the nonrunning coupling X(x) = A*, we have the an- 
alytical result £ = — l/(87r) = —0.0398 corresponding 
to Eq.(g]), P = -2o) 3 /(3Cf), which agrees with the nu- 
merical result £ = —0.0400 in the vicinity of the critical 
coupling. Incidentally, in the case of the two-loop cou- 
pling we may take a fitting function other than the above: 
£(1/A cr - 1/A*) 3 / 2 with £ = -0.0155, which yields much 
better fitting. 

How about the exponent of the nonperturbative beta 
function with respect to Col Let us write 



la {Gl 



/(A*)£ K , 



ID 



\1V I 



where /(A*) represents some function of A* and the in- 
stantaneous exponent k = k(A„) is extracted from the 
relation 



~ a 



OJ 



( 4 «) 



(82) 



We depict the numerical results in Fig. [13J Near the 
critical edge, the numerical values of n = k(A*) are 



k~ 2.73, 2.75, 2.95, 



(83) 
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FIG. 14: Conjecture of the shape of the /3 function including 
the nonperturbative behavior. The bold solid and dashed 
curves correspond to the conjectured j3 function and the 
purely perturbative one, respectively. 

for the two-loop exact solution, the parabolic approxima- 
tion and the fixed coupling case, respectively. Because 
the linear and multiple zeros of the beta function with 
respect to a(/i) = a ~ a* correspond to k|a»=a„ = 2 
and k|a„=a ci: > 2, respectively, the numerical results ob- 
viously show that the nonperturbative beta function has 
the multiple zero at the critical edge. If we smoothly ex- 
trapolate the behavior of k to the criticality, the behavior 
of the nonperturbative beta function at the criticality will 
be (3 cx — w 3 , in accord with the above least-square fitting 

It is noticeable that the nonperturbative beta function 
has a multiple zero at the critical coupling a ~ a» = a cr , 
as shown in Eq. (|80[) . which corresponds to the essen- 
tial singularity scaling Eq. (|58|) . On the other hand, it 
is not the case in the perturbative (two-loop) expression 
([2"12]) . which has a linear zero, (3 ~ a — a* ~ a cr — a* at 
criticality a = a cr . Therefore the actual beta function 
is crucially different from the perturbative beta function 
which should be modified in the IR region where the non- 
perturbative dynamics responsible for the mass genera- 
tion is dominant. The full j3 function including the per- 
turbative and the nonperturbative contributions is thus 
suggested in Fig. Q3J We hope that the lattice studies 
will clarify this nature. 

Now we discuss the behavior of the TD mass in the 
limit toward the criticality. Through the PCDC, Eq.©, 
the vacuum energy is connected with the TD mass, 

FtbMtb = -de(0£) = -AdoV, (84) 
where -Ftd, Mtd and dg(= 4) represent the TD de- 



cay constant, the TD mass and the scale dimension of 
the trace of the energy-momentum tensor. Our results 
qualitatively agree with the conclusion in the case of SD 
equation with the nonrunning gauge cou plin g, Eq . ([T12]) . 
which is in disagreement with Refs. (2I |27|: There is 
no true (massless) NG boson for the conformal sym- 
metry at the criticality, unless the TD decay constant 
diverges, i.e., gets decoupled. Such a decoupled TD 
was in fact implied by the idealized limit of the holo- 
graphic TD [23. In the realistic situation of the TC 
model building m/Axc is not arbitrarily small but only 
m/Axc ~ to/Aetc ~ 10~ 3 — 10~ 4 , so that the "mass- 
lessness" and "decoupling" are somewhat milder. In the 
quantitative sense our results, though valid only in the 
vicinity of criticality, indicate rough idea about the mass 
and decay constant of TD as follows. Substituting our 
numerical result, Eq. (|77)l. into Eq. (1541 . we find 

4 

A/| D ~ 3.02 77^-. (85) 

^TD 

Furthermore, by using F^/m ~ 0.41 */ NjjNtc /3 in 
Eq. ([7I]). we obtain 

Fir Ftd \ V 2N D V 2N D V N TC J 

The TD with mass, say Mtd ~ 500 GeV, would require 
TD coupling smaller than that of the Standard Model 
Higgs by -FV/Ftd — 3/5 up to model-dependent factors 
of Nj,Ne> and A^tc besides other dynamical details. If 
the TD mass is much smaller, Ftd <§C F w , on the other 
hand, it could lead to a decoupled TD, which might be a 
candidate of dark matter. Detailed studies are required 
in order to confirm whether or not such a decoupled TD 
satisfies conditions for dark matter. 



IV. SUMMARY AND DISCUSSIONS 

We have studied analytically the improved ladder SD 
equation with the parabolic approximation for the beta 
function and also analyzed numerically the solution of 
the ladder SD equation with the two-loop exact running 
gauge coupling. 

We explicitly calculated the techni-gluon condensate 
near the conformal edge and found that the behavior 
is (G 2 ill )/m A ~ (a/a CI — 1)~ 3 / 2 — > 00 (a — ¥ a cr ) with 
m <C Atc, in accord with Ref. [23[ but in disagree- 
ment with the assumption of Refs. 0, H3- The nu- 
merical calculation is consistent with this analytic re- 
sult. This situation is different from QCD having no 
approximate conformal symmetry, where Aqcd ~ m and 

(Cpi/)full ~ (G^J) perturbative ~ Aqq D . 
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On the other hand, the vacuum energy (divided by to 4 ) 
is finite, V "/to 4 — > const, even in the critical limit, as in 
the case of the fixed gauge coupling. Our result for the 
vacuum energy only yields a combination of the mass 
Mtd and the decay constant Ftd through PCDC but 
not each of them separately, as was the case for most of 
the discussions in the literature. Combining the PCDC 
relation, Eq. ([84]), F^ D M| D = -4(0£), with the numer- 
ical result for the vacuum energy, Eq. (|77p . we found 
M 2 D = 3.02?7to 4 /F 2 d with Tj = N TC N f /(2n 2 ), Ea.([85]t. 
This relation implies Mtd/w ~ m/i*TD near the confor- 
mal edge and hence naturally Mtd = 0(m) in contrast 
to Refs. [27ll28|. (A similar conclusion was made in a dif- 
ferent context |38|.) As an extreme possibility we could 
have Mtd/w — > only when to/Ftd — > and the TD 
gets decoupled. If such an idealized decoupled massless 
TD is realized at the conformal edge, the light decoupled 
TD as a pseudo NG boson slightly off the conformal edge 
may be a candidate for the dark matter. 

The scale anomaly formally yields the relation 
/3/(4a) = (00 1 (G 2 „), so that the above results imply the 
nonpcrturbativc behavior of the beta function reflecting 
the nonpcrturbativc effects of the dynamical mass gen- 
eration, f3(a) ~ — (a/a cr — l) 3 / 2 . Numerically, we ob- 
tained /3/(4a) = £(a/a cr - l) 3/2 /A* with £ = -0.026 
for the two-loop running gauge coupling. The absolute 
value of the coefficient is smaller than that for the fixed 
gauge coupling, £ = —0.04. However, the exponent k of 
the nonperturbative beta function at the conformal edge 
seems universal, n = 3, where /3/(4a) oc (a/a cr — 1) K / 2 . 
This nature of the nonperturbative beta function hav- 
ing the multiple zero is crucial to reproduce the essential 
singularity scaling Eq. (|58|) . 

We have settled some of the controversy related with 
the TD mass raised within the improved-ladder SD equa- 
tion. However, several issues remain to be explored: 

In particular, a central problem is how large the TD 
mass Mtd is. In order to discuss collider phenomenology 
of the TD and also check whether or not a decoupled 
TD is in fact realized near the conformal edge, we should 
obtain mass Mtd and the decay constant Ftd separately. 

Thus we would need more information other than the 
vacuum energy. As was mentioned in the Introduction, 
such a calculation was in fact done in a most straightfor- 
ward way [2l| , based on the SD and BS equations in the 
improved ladder approximation with the two-loop run- 
ning coupling constant having the CBZ IR fixed point, 
which suggests Mtd ~ V2m, Eq. (fl"4|) . without evidence 
of the decoupled light TD much smaller than to. Note 
however that this calculation was actually done only nu- 
merically and at slightly off the conformal edge, and 
hence the result is not conclusive about the very close 
point to the conformal edge. 

On the other hand, in the holographic framework [23[ 



which has a wider parameter space than that of the (im- 
proved) ladder approximation so as to adjust the S pa- 
rameter arbitrarily small, it was shown that at the limit 
of conformal edge to/Atc — * the techni-gluon con- 
densate vanishes r — > oo, with T parameterized as in 
Eq. ([17[) , which in turn implies Mtd/™ — ^ at the sacri- 
fice of decoupling to/^td 0, although such an extreme 
case is unlikely for the realistic setting of the typical TC 
model building (slightly away from the conformal edge 
to/Atc = 10~ 3 — 10 -4 ) where the holography suggests 
M TD /m = 0(1), or M TD ~ 600 GeV (Eq.®). So al- 
though the theoretical possibility for the decoupled TD 
at the conformal edge is not completely excluded, there 
is no signature of such a possibility at least in near con- 
formal edge region relevant to the realistic TC model 
building. 

We have not included interactions like ETC (Ex- 
tended Technicolor) between the tcchni-fcrmions and the 
quarks /leptons which should be included to give mass to 
the quarks and leptons in the realistic TC models. In- 
cluding these interactions also induce additional interac- 
tions among the techni-fermions themselves, which may 
be described by the effective four-fermion interactions in 
addition to the TC gauge interactions we have discussed 
("gauged NJL model"). Such ETC-like effects on the TD 
mass were already studied intensively in the ladder SD 
equation with nonrunning gauge coupling [l(| EH Il3j |. 
with the results Mtd ~ C( m ), i.e., against very light 
TD mass and decoupled TD, as was mentioned in the 
Introduction. As is clear from our arguments in this pa- 
per, the situation with the additional four-fermion in- 
teraction in the ladder approximation with nonrunning 
gauge coupling will remain essentially the same in the 
improved ladder SD equation with the two-loop running 
gauge coupling. Moreover, more elaborate calculations 
in the gauged NJL model [l5[ showed that Mtd —> V^rn 
(Eq.([T2])) in such a limit along the whole critical line 
(0 < a < a CI ). Note also that the result Ea-tfTH)) QJ] 
is consistent with Eq. flT4l [2lj which is the straightfor- 
ward computation of the spectra within the ladder SD 
and BS equations in the improved ladder approxima- 
tion with the two-loop gauge coupling near the conformal 
edge (without four-fermion coupling) . Note however that 
these calculations did only for the inverse propagator at 
zero-momentum but not the pole mass (on-shell mass), 
and hence are still not conclusive. We definitely need 
more reliable calculations such as the lattice simulations 
about the TD spectrum. 

We have considered TD as a bound state of techni- 
fcrmion and anti-techni-fcrmion both of which acquired 
mass to. The mass to spontaneously breaks chiral sym- 
metry and at the same time breaks spontaneously and 
explicitly the scale symmetry, the scale anomaly due to 
this mass generation being of order C(to 4 ) as we com- 
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puted explicitly in this paper. Hence such a bound state 
should have mass of order C(m)(<C Atc)- O n the con- 
trary, it was argued [23[ that the techni-glueball mass 
should be of order O(Atc), since the scale-symmetry 
breaking free from the techni-fcrmion mass generation 
is due to the scale anomaly of order 0(A TC ) correspond- 
ing to the usual perturbative running of the coupling for 
fi > A T c ({00 perturbative in Eq.((35]) ). Then we expect 
little mixing between our TD and the techni-glueball, in 
sharp contrast to QCD where the flavor-single scalar me- 
son (analogue of TD) and the scalar glucball may mix 
strongly. More reliable calculations are of course highly 
desired. 

In this paper, we assumed that the fermion loop is 
dominant in the techni-gluon condensate and the vac- 
uum energy. In principle, there might exist nonpertur- 
bative techni-gluonic effects. It is difficult to estimate 
it in our approach, however. A lattice simulation may 
also resolve this issue. We shed a light on the problem 
which has made confusion in the improved ladder approx- 
imation with the two-loop running gauge coupling. We 
clarified it in the analytical and numerical ways within 
the same framework of the improved ladder SD approx- 
imation. There should exist nonperturbative effects be- 



yond the improved ladder approximation, although it is 
very unclear whether or not they are relevant. Toward a 
conclusive answer, one might challenge a more rigorous 
approach such as a lattice gauge theory. We hope that 
our results will be reconfirmed more rigorously in future. 

Needless to say, our analysis in this paper is also ap- 
plicable to dynamical symmetry breaking scenarios with 
large anomalous dimension/ conformality other than TC, 
such as the top-mode standard model with extra dimen- 
sions which has a UV fixed point [36[ , higher representa- 
tion quark condensate model \^ , or even the QCD with 
finite temperature where the running of the coupling will 
be frozen (mocking conformal) in the IR region below the 
temperature scale, etc.. 
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